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(57) Abstract 

A method and apparatus for detecting amplitude variations in the rhythmic oscillations of a physiological response 
pattern in a frequency range of interest A sensor (I) is used to detect an occurrence of an event in the cycle of a physiolog- 
ical response and the interval between each reoccurring event is timed and placed in a buffer (9). The output from the buf- 
fer (9) is read out at predetermined time intervals and fed to a filter (90) which determines the aperiodic portion of the sig- 
nal and substracts that aperiodic portion to output a residual data signal. The residual data signal is fed to a band pass 
filter (94) which filters in a region determined by the predetermined frequency range of interest and which outputs to a cal- 
culation and display device (96) wherein the variance of the rhythmic oscillation is calculated and displayed. 
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Description 

Method And Apparatus For Evaluating 
Rhythmic Oscillations In Aperiodic 
Physiological Response Systems 

Technical Field 

The present invention relates to the 
quantification of the variability of physiological 
activity and particularly the quantification of 
oscillations in physiological response patterns such as 
the heart rate. This invention was made with 
Government support under National Institutes of Health 
grants Nos. KO2-MH-0054; MH-18909; HD-15968; and HD- 
05951. The Government has certain rights xri this 
invention . 

Background Art 

The measurement and determination of effective 
diagnostic information from the output of physiological 
response systems (e.g., heart period patterns/ 
peri phera l vasomotor activity/ electrodermal 
potentials/ electric potentials from the scalp such as 
EEGs / blood pressure/ temperature and all other 
physiological activity which may be indexed by time) is 
complicated by the nature of these response patterns 
which often are characterized by rhythmic oscillations 
superimposed on an aperiodic base line. Extremely 
complex underlying mechanisms of human physiology 
underlie the mechanisms which produce the signals. 
That is/ physiological response systems tend not to be 
determined by a single input but are the result of 
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complex interactions of numerous, often undefined/ 
mechanism. For example, the nervous system has a 
profound impact on many physiological responses by 
modifying "horaeostatic" oscillations which represent, 
in particular circumstances/ known physiological 
mechanisms. More specifically, the heart period (the 
time between successive heart beats) presents 
oscillations which are located in frequencies common to 
other physiological response systems. The heart period 
oscillates at the breathing frequency and at the 
frequency at which blood pressure and peripheral 
vasomotor activity also oscillate. The oscillations in 
the heart period at the "respiratory" frequency and at 
the "vasomotor" frequency may be interpreted as a 
15 indicaricn of specific physiological mechanisms. 

When studying the measured physiological activity 
in terms of oscillations, the parameters of interest in 
order to provide information of diagnostic value, 
include the amplitude of the oscillation, the phase of 
20 the oscillation relative to other periodic 

physiological functions at the same frequency and the 
coupling or coherence between two or more physiological 
systems at the same frequency. There are numerous 
methods used to separate the signals of interest or to 
25 de trend physioi-ogical data. Many of these on-line 
devices for monitoring the physiological response 
system s inclu dehigh-pass, low-pass or bandpass 
^£i-Ltera. other methods— include specific statistical 
analysis which have been developed for engineering and 
economics applications but which are attempted to be 
used in regard to the detrending of physiological 
data. in general, most of these prior art procedures 
assumed that the trend which is being removed may be 
characterized by a linear regression or the sum of slow 
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sine waves. While these methods appear to function 
well in some areas of physiological monitoring 
including respiration and the electrocardiogram/ they 
are rather limited to those instances where the 
variance associated with the rhythmic oscillations 
being studied is large relative to the instability of 
che base line upo n which these oscil lations are super- 
iraposed^ — /On theother hand/ when the variance 
asociated with the oscillations of interest is 
10 extremely small relative to the total variance of the 
physiological response system/ then the above 
assumptions associated with the previous methods and 
apparatus of detrending do not apply. 



An example of an instance where the prior art type- 
15 of filtering for purposes of detrending the 

physiological data fails is that of the amplitude of 
fetal heart rate oscillations which are very small 
relative to large changes in the heart rate which have 
been associated with uterine contractions. This is 
20 especially pertinent when a compromised hypoxic fetus 
exhibits massive heart "rate shifts in the base line in 
response to the uterine contractions/ thus making it 
very difficult to accurately estimate the fast periodic 
heart rate activity which is of diagnostic value. 

25 Most of the statistical procedures which are used 

to assess the characteristics of periodic processes 
such as the amplitude of rtrytHmlc oscl 1 h-tlons lnvolve 
attempts to detrend the base line to remove aperiodic 
components from the data set* These periodic processes 

30 which are embedded in the complex signal are attempted 
to be removed through the use of sequence involving 
detrending, filtering, and describing the amplitude and 
periodicity with spectral analysis such as the fast 
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fourier transform ( FFT) • The detrending and filtering 
produces a "processed" signal by removing the aperiodic 
component which allows for the use of statistical 
procedures to evaOuajce^the^ tAe rhythmic 



J3S cili a ti on s^ -\Thi3 -processed" signal is decomposed 
through the use of spectral analysis and the variance 
is partitioned into constituent frequencies. That is, 
the variance is described as the sum of sine waves of 
various amplitude and frequencies. The problem with 
this process is that it may result in faulty 
interpretations of data if the data set being processed 
violates specific statistical assumptions necessary for 
proper interpretation. 
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Spectral analysis may be used to accurately 
identify and quantify periodic components in 
physiological response systems when they are only 
minute base line shifts or when the base line trend can 
be easily removed prior to analysis. Spectral analysis 
assumes that the data set being analyzed is weakly 
stationary. a process is weakly stationary, if its 
mean and variance are independent of time and its 
autocovariance function depends only on lag (c. 
Chatfield, The Anal ysis of Time Series; Theory and. 
Practice,, .Chapman and. Hall 197-5-)-. spectral analysis 
provides reliable and interpretable estimators of the 
amplitude of a periodic oscillation only if the data 
are at least weakly stationary. 



Another of the unfortunate physiological response 
system characteristics is that they are not 
30 "stationary". This means that physiological response 

systems are not perfectly sinusoidal and that they have 
complex shifts in both the mean level and the 
variance. Thus, by their nature, they violate the 
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assumption of stationarity. Quite obviously then, the 
spectral analysis to evaluate the amplitude of rhythmic 
oscillations will result in unreliable estimates of the 
amplitudes of the rhythmic process at any specific 
frequency band. By appropriately removing the complex 
base line trend/ it would be possible for the amplitude 
of the periodic oscillations to be accurately measured 
in the filtered data set; however/ all of the existing 

filter! ng met hods and_d ejsaces.-w.hxch-have--been M uaed.t.o 

"de trend" physiol ogical response activity in order to 
remove the" shiftin g base line have ,made_f Aulty 
assumptions, _ ; g 

Many existing physiological monitoring devices 
such as polygraphs/ electroencepniographs, and 
electromyography have hardware filters which function 
as high— pass / low-pass / or ban dpass filters . As 
previously discussed/ this reflects an assumption that 
there are no aperiodic components in the data and 
merely that the filters pass the frequency band of 
interest to the output. Unfortunately/ since the data 
of most physiological systems contain aperiodic trends / 
the amplitude of the frequencies passed by the various 
flljters wiJLl be partially a function of the amount of 
variance passed through the filters which la, In 
reality/ a portion of the complex aperiodic trend 
discussed above. Thus/ in essence/ the hardware 
fil.te-rs of... prior art devices assume that the base lines 
are merely the sum of slow sine waves and a linear 
trend. If the trend is complex and cannot be described 
by a linear trend or a sura of sine waves with known 
periodicities and amplitudes, then the sine waves 
necessary to describe the slow complex trend may 
include faster periodicities superimposed on the 
frequencies of interest. Therefore/ the operators of 
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the device must know beforehand the shape of the trend 
to be subtracted from the data set. In the case of 
spectral analysis, it would be necessary to subtract 
the spectral densities associated with the trend from 
the spectral densities associated with the total data 
set. This is totally impractical because, with the 
filters being used, the operator would never be able to 
separate what component of the variance being passed by 
the filter is associated with the trend from that 
component of the variance which would be associated 
with the periodic process. Moreover, it would preclude 
the ability of the operator to monitor the changing 
conditions of the periodic physiological process in an 
on-line operation. 

Other methods of operating a filter for removing 
trends include the use of what is called "successive 
differencing^. This method consists of successively 
subtracting values through the entire data set 
involving, for. example, the subtraction of data point 
number 1 from data point number 2 and data number 2 
from data point 3, etc. Due to the transfer function 
of this filter the method may result in an 
underestimate or overestimate of the spectral densities 
depending on which frequencies are of interest to the 
investigator and thus may result in a contamination of 
the estimates of the variance at any specific 
frequency. Moroever, the "successive difference" 
filter is similar to linear detrending and- suffers from 
the same problem of passing' variances in higher 
frequency trends which are components of the aperiodic 
trend. Other methods include low order polynomial 
detrending techniques which do not succeed in removing 
the trend and which also result in an alteration in the 
shape of the spectrum by influencing both the 
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Identification of the peak frequency and the estimate 
of amplitude at a given frequency. 

The clinical and diagnostic value of overcoming 
the errors brought about by incorrect or simplified 
assumptions in the prior art devices discussed above 
can be particularly seen in a specific situation where 
the amplitude of the oscillation of a physiological 
process may serve as an indexing variable of a specific 
underlying mechanism. For example, in the case of a 
heart period , it is possible to interpret the amplitude 
of the oscillation of the heart period at the 
respiratory frequencies or respiratory sinus arrhythmia 
( RSA) as an index of the influence of the vagus (10th 
crainial nerve) on the heart* Briefly stated/ the 
r espr ra: t o ry system transmits afferent information to 
the brainstem where it "gates" (turns off and on) the 
vaga~L efferents to the heart (i.e./ vagal efferent 
activity is reduced during inhalation and reinstated 
during exhalation). Thus/ with regard to heart period 
oscillations occurring in the respiratory frequencies/ 
the amplitude of such* heart period oscillations conveys 
information regarding the "vagal tone" effect on the 
heart. With regard to clinical and diagnostic 
relevance/ higher order .central nervous system 
disorders . such as intracrainial hermorrhage result in a 
decrease in the vagal efferent influence on the 
heart. Therefore/ the amplitude of the heart period 
oscillation in the respiratory frequency band (RSA) may 
provide a "window to the brain" and an early assessment 
screening of the central nervous system dysfunction. 

One of the more important points concerning the 
relationship exhibited by a fetal heart period is that 
the heart period pattern exhibits small oscillations at 




WO 84/02641 



PCT/US83/02031 



-8- 

the periodicities associated with breaching- in the 
newborn which frequencies are most likely 
representative of the RSA in the fetus. Because the 
periodicities account for much less than one percent of 
5 the total variance of the heart period pattern, the 

above-discussed points with regard to prior art methods 
of evaluating and detecting the amplitude of the 
periodic function become more critical and more prone 
to error because the percent of variance that the 

0 specific periodic function accounts for is extremely 
small relative to the total variance of the 
physiological response pattern in the fetal heart 
period. 

In the fetus, the heart period- is mainly 
5 influenced by the f e to-maternal movement and the impact 
of. uterine contractions during labor. Spectral 
analysis of fetal heart period utilizing any one of the 
above-discussed filtering techniques would mask the 
presence of the small oscillations because the percent 
3 of variance would be so low that it would not result in 
a significant or even recognizable peak in the spectrum 
and because the variance from the complex trends 
related to either f eto-ma ternal movement or uterine 
contractions would produce more variance in the 

1 frequency band of RSA than the RSA itself. 

In summation then, the spectral analysis provides 
interpretable estimates of the variance (amplitude) on 
specific frequency bands only when the data do not 
violate specific assumptions. Most data sets derived 
) from physiological systems such as heart period 

activity contains aperiodic components and violate the 
critical assumption of weak stationarity (i.e., the 
expected values of the mean and the amplitude are 
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independent of time) • This is critical when the 
variance of the oscillations of periodic physiological 
activity conveys critical information regarding the 
condition of the organism as the instance when the 
5 reduction in the amplitude of RSA in the human neonate 
is associated with brain damage and/or nervous system 
conditions threatening the viability of the infant. 
Likewise/ spectral analysis and most of the filtering 
techniques which attempt to remove the trends/ are not 
10 readily adapted for rapid on-line use. 

Disclosure of the Invention 

It is a primary object of the present invention to 
provide a device which incorporates a filtering system 
which "on-line" removes the variance associated with 
15 complex trends and slow sine waves enabling accurate 

estimates of more rapid periodic oscillations when they 
are superimposed on complex trends and slow sinusoids* 

It is also an object of the present invention to 
incorporate a procedure which provides an accurate 

20 evaluation of the amplitude of a rhythmic oscillation 
on a frequency band when the rhythmic oscillation is 
s upe r imposed "on an aperiodic response pattern - . The 
procedure involves the removal of the effect of base 
line drift by detrendl ng the ba se line activity with a 

25 "moving polynomial filter* (MPF) . The moving 

polynomial filter consists of two stages: first, it 
smooths the base line pattern by fitting a piece-wise 
polynomial to the base line; second the smoothed base 
line pattern is subtracted from the original data 

30 set. The residuals represent a filtered data set with 
a low frequency cut-off. The output of this procedure 
is then bandpassed to allow only the variance of the 
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data set associated with che frequency of interest to 
pass. The output of these two filtering devices and 
methods provides the component of variance associated 
with rhythmic oscillation of interest. 

It is also an object of the present invention to 
provide a moving polynomial filter (MPF) which provides 
an ability to assess on-line the continuous shifting 
amplitude of the oscillation of interest during 
situations when the periodic process being monitored is 
not weakly stationary and which functions to assess the 
amplitude of oscillations on-line without the 
collection of data for an entire session prior to 
detrending. 

It is also an object of the present invention to 
provide an apparatus and method for detecting central 
nervous system dysfunction in humans by providing an 
arrangement of detectors and filters which output a 
clinically usuable signal providing an indication of 
disorders in the central nervous system through the 
monitoring of the heart rate. 

Brief Description of the Drawings 

A more complete appeciation of the invention and 
many of the attendant advantages thereof will be 
readily obtained as the same becomes better understood 
by reference to the following detailed description when 
considered in connection with the accompanying drawings 



FIGURE la is a representation of a physiological 
response pattern with a rhythmic oscillation of an a 
priori frequency band superimposed on an aperiodic 
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response pattern. 

FIGURE lb is a representation of the smoothed base 
line template to be removed from FIGURE la. 

FIGURE lc is the residual from resulting 
subtracting of the template of FIGURE lb from the 
pattern FIGURE la. 

FIGURE 2 is a block diagram of the device for 
evaluating the amplitude of the rhythmic oscillations 
superimposed on the aperiodic response patterns. 

FIGURE 3 is a graphic illustration of the moving 
polynomial filter of FIGURE 2. 

FIGURE 4 is an illustration of the clinical 
importance of the method and apparatus of the present 
invention of FIGURE 2. 

Best Mode for Carrying Out the Invention 

Referring now to the drawings/ wherein like 
reference numerals designate identical or corresponding 
parts throughout: the several views/ and more 
particularly to FIGURE la thereof/ there is shown a 
superimposition of the amplitude of a rhythmic 
oscillation of a known frequency band on a aperiodic 
response pattern. This com pi-ex-a±g.na]=sseEveav.«-a 
model for physiological response patterns. 

A specific example of the type of complex system 
illustrated in FIGURE la is that associated with the 
oscillation of heart period which must be "separated" 
in order to provide the measurement of the vagal tone 




WO 84/02641 



PCT/US83/02031 



-12- 

which, as previously discussed, is the spectral 
representation of the amplitude of periodic" heart 
period activity associated with the "gating" by 
respiratory activity of the vagal efferents to the 
5 heart- The accumulation of spectral density estimates 
of heart period activity associated with the 
respiratory frequency band provides an accurate 
measurement of respiratory sinus arrhythmia (RSA) . The 
vagal coney as measured by the amplitude of the RSA, 
10 provides a key to the link between the heart rate 

response and an underlying physiological dysfunction* 
More particularly , the amplitude of the RSA may be 
employed as an index of central nervous system 
dysfunction or influence* 

15 In order to provide clinically reliable 

information for use as a diagnostic tool, the rhythmic 
oscillation of FIGURE la must be separated from the 
underlying aperiodic signal. 

Thus, if the heart period activity is filtered to 
pass only the frequencies associated with respiration, 
the new heart period series would represent the RSA. 
Because the amplitude of any sinusoid is mono tonically 
related to its- variance, the amplitude of the sinusoid 
resulting from the bandpass filtered heart period 
activity would represent a sort of breath-by-breath 
RSA. An instantaneous estimate of the vagal tone could 
be derived by linear interpolation connecting the peaks 
of the adjacent filtered heart period oscillations. 

In physiological response patters, non-sinusoidal 
30 proceses distort the estimates of the underlying 

rhythmic processes. Thus, it is necessary to remove 
these aperiodic influences prior to estimating the 
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variance or amplitude associated with specific rhythmic 
processes. Generally/ the methods used to detrend the 
base line drift have assumed that aperiodic influences 
may be removed by detrending with a linear regression 
5 or fitting the entire data set to a lower-order 
polynomial. However/ most base line trends in 
physiological response patterns are aperiodic so that 
neither a polynomial fit not a sum of slow sine waves 
can adequately represent the trend. For example/ in 

10 the human fetus / there is a relatively fast frequency 
(25 to 75 cycles per minute) in tne fluctuation of the 
heart period pattern. This pattern is difficult to 
evaluate because it may be superimposed on slower 
shifts in the heart period. One source of the slow 

15 shifts is the mother's uterine contractions. Thus/ the 
rhythmicity of the fetal heart p erio d pattern/ which at 
specific times is stationary/ is superimposed on a slow 
complex trend . 

The complex trend in the fetal heart period 
20 pattern as related to the uterine contractions cannot 
be removed by linear detrending nor by detrending low- 
order polynomials. However/ a local cubic polynomial 
moving average may be stepped through the data to 
produce the desired- series. The moving average- is a 
25 weighted sum applied locally along the curve to obtain 
a smooth curve or template as shown in FIGURE lb. This 
template of slow activity is then subtracted from the 
raw data— and— th^ residual s^ri:es of FIGURE lc contains 
the rhythmic process which is free from the influences 
30 of the aperiodic component and is stationary over short 
segments. In the case of the fetal heart period 
pattern/ the residuals are a time series consisting of 
the fast fluctuations. These fast fluctuations may not 
be stationary over long periods of time/ because the 
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amplitude of the oscillations may change as. the 
condition of the fetus's central nervous system 
changes - 

The FIGURE 2a, 2b details an apparatus which 
functions to remove the rhythmic oscillation of 
interest from a physiological response pattern as shown 
in FIGURE la, lb, lc. The sensors 1 are provided by 
pickups which sense physiological activity from che 
body to include electrodes to detect endogeneous 
bioelectric potential as the beating of the heart in an 
EKG as well as the volumemetric and plethysmography 
sensors to assess response systems such as respiration 
and blood flow. The output from the sensors may be 
either seen to be a point process (FIGURE 2a) which is 
15 manifested as events in time such as the beating of the 
heart or a continuous process (FIGURE 2b) which is 
manifested as continuous changes over time such as the 
changing the circumference of the chest during 
respiration and the electrical potential of the finger 
20 or scalp. In using a measurement of a point process as 
shown in FIGURE 2a , the component of the sensor device 
which is called the "event or peak detector 1 * 5 detects 
the. time of occurrence of specific events in a point 
process. In the instance of heart perio_ds this could 
25 be the detecting of an occurrence of a R wave in the 
EKG representing the contraction of the ventricles of 
the heart. The output of the peak detector is fed to a 
clock 7 which times the irttBCTals- between successive, 
events. 

30 The values of the clock are stored in a buffer 9 

for later use and the output of this buffer 9 is fed to 
an outlier compensator 11 which monitors the timing 
between events of interest such as the occurrence of 
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the R waves. The compensator contains information 
concerning the expected range of time intervals between 
the R waves based upon a knowledge of the subject being 
tested. This outlier compensator maintains a specific 
5 interval between events occurrences which are later to 
be analyzed by the remainder of the circuitry. 

Some of the problem overcome by the outlier 
compensator include the problem which occurs with 
regard to specific ventricular arrhythmias which result 

10 in an out-of-place event detection. For example, 
ventricular extrasystoles are followed by a 
compensatory pause. Ventricular extrasystoles do not 
interrupt the regular discharge of the sinoatrial 
node. me duration of the beat with the ventricular 

15 extrasystole plus the beat following with the 

compensatory pause are equivalent to the duration of 2 
normal beats. Of importance is the fact that centrally 
mediated rhythms / such as RSA, which are manifested in 
the heart period pattern are naturally occurring 

20 arrhythmias of the sinoatrial node. Ventricular 
extrasystoles do not interrupt or influence these 
rhythms. Thus / to provide ^an accurate estiraate^of the 
variances associated with heart period rhythms such as 
RSA~> it is-nccassar-y t©~ .correct the ventricular 

25 extrasystoles and compensatory pauses. The outlier 
compensator would adjust for the ventricular 
extrasystole and compensatory pause by adding the two 
time periods together and dividing by two. Thus the 
outlier compensator is no more than an adder and a 

30 divider circuit arrangement programmed to operate when 
the duration between events fed from the buffer 9 is 
outside of a range of expected event intervals based 
upon the subject being tested. In this isntance/ the 
outlier compensator removes this event by adjusting the 
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space between the events in accordance with its 
expected spacing intervals. Other areas where the 
outlier compensator 11 is functionally important 
include the instance where there is a missing beat due 
5 to the failure of the components of the "event 

detector* to be triggered or in the instance where it 
is triggered prematurely resulting in too many beats. 
In each of these instances the outlier compensator 
checks and maintains a proper predetermined interval 
10 range between the beats by either adding together the 
intervals or dividing an interval by two. 

The output from the compensator is fed to the time 
window 20. This time window 20 measures the time 
interval between the event being detected such as the R 
15 wave within a specific timeframe or window. For 

example, if a time window operates every 500 msec then 
the data set is transformed into an estimation of heart 
period for each 500 msec time window and is outputted 
to the moving polynomial filter 90. The determination 
20 of the time interval is in cooperation with the 

frequency characteristics of the underlying process and 
the moving polynomial filter 90. The timing or the 
duration of the time window is critical to prevent 
aliasing. As^is known in" the art of signal ^ceasirig 
25 aliasing occurs when there are variances associated 
with frequencies which are faster than those 
frequencies being detected and those variances from 
outsride -the^r^qTie^y-bami-of interests are folded back 
(i.e., added to) onto those frequencies of interest. 
30 In order to protect against this aliasing in the heart 
rate example the time window is set to be such that the 
time window samples are approximately twice as fast as 
the R wave events which identify each heart beat. Thus 
if the normal heart rate is assumed to be 60 beats per 
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minute the sample or time window is set at twice the 



window operates every 500 msec to dump its data output 
to the filter. The use of a sampling rate set at twice 
the frequency of interest is based upon a statistical 
rule set down by Blackman and Tukey (R.B. Blackman & 
J.W. Tukey/ The Measurement of Power Spectra / 1959/ New 
York: Dover) which basically realizes that if the 
sampling is twice as fast as the event in being . 
detected it is sufficient to provide all of the details 
of the event without any aliasing from frequencies not 
of interest. In other words it prevents the folding 
back phenomena. Quite obviously the time window could 
be set to smaller intervals (i.e./ faster sampling 
rate) as long as it obeys the basic rule of sampling at 
least twice the frequency of the event of interest. 

To adequately employ a moving polynomial filter 
(MPF)/ the data must be sampled at appropriate time 
intervals. Since the moving polynomial filer (MPF) 
functions in the time domain and assumes that the data 
are sampled at equal time intervals/ it is important 
that those of the processes which are called 
"continuous processes" must be sampled at a rate to 
generate discrete processes containing the variance 
associated with the fastest oscillations in the 
process. The choice of the sampling interval is 
critical because it is not possible to identify rhythms 
whose period is less than twice the time between 
observations. It can bs mathematically demonstrated 
that faster frequencies (i.e./ faster than twice the 
duration of the sampling interval) will be "folded 
back" or "aliased" on the slower frequencies. This may 
result in an overestimate of variances in the frequency 
band of interest (P. Bloomfield/ Fourier Analysis of 



speed or 2 events per second. In other words, the 
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Time Series: An Introduction , Wiley/ 1976) , The 
FIGURE 2b details the apparatus necessary for a 
continuous process wherein the same numbers represent 
identical apparatus in FIGURE 2a. Basically/ the 
5 Analog to Digital convertor 14 is used instead of the 
Peak detector 5/ the clock 7, the Buffer 9/ the 
compensator 11 and the time window 20 due to the above 
discussed nature of a continuous process. The sampling 
rate of the A/D convertor 14 is analogous to the' time 

10 window 20 of the Point Process because of the need to 
sample a "Continuous Process" at a rate so that the 
variance associated with the faster oscillations are 
contained in the "discrete processes" from the A/D 
convertor 14. It also should be noted that continuous 

15 processes must be sampled at equal time intervals while 
point processes describe by interevent time intervals 
(e.g./ heart period: the time between successive 
heartbeats) must be transformed by appropriate 
weighting or sampling into equal time intervals. The 

20 prior art does not provide for filtering point 

processes in the time domain. Accordingly, the device 
of the present invention tranaforras the point process 
to a discrete process sampled at equal time 
intervals* This enables the application to time domain 

25 filters to point processes. With regard to continuous 
processes/ it is necessary to have prior knowledge of 
the physiological system being evaluated. If it is not 
known what the periodic components embedded in the 
physiological process are, aliasing may result in the 

3 0 estimates of variance being uninterpre table. The 

moving polynomial filter (MPF) provides furthermore a 
facilitation in the ability to assess the amplitude of 
oscillations on-line because it does not necessitate 
the collection of data for the entire session prior to 

35. detrending. The MPF also functions in situtions in 
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which Che rhythmic process being studied is not 
stationary and provides the ability to assess on-line 
the continuous shifting amplitude of the oscillation of 
interest. 

5 The moving polynomial filter functions as a high- 

pass filter when the smooth template is subtracted from 
the data set. The frequencies passed by the filter are 
a function of the duration of the polynomial which is 
defined by the number of data points incorporated into 
10 the moving polynomial with each data point representing 
a time— sampled measure. In physiological applications 
a third order polynomial is most appropriate because it 
enables the removal of cubic dips in the trend. 

The moving polynomial filter 90 which is fed with 
15 the time windowed data set consists of a digital filter 
which is stepped through the data set. A low order 
polynomial is fit on a "local" level to describe the 
aperiodic trend. Although a polynomial fit or trend 
over the entire data set suffers from many practical 
20 and statistical problems/ any smooth function can/ 

under very general conditions/ be represented locally 
by a polynomial to a fairly high degree of accuracy. 
By stepping a localized polynomial through the data/ it 
is possible to smooth the data set and- to describe the 
25 complex aperiodic trend. When the lt stDooched' , trend is 
subtracted from the original data/ the residuals 
represent a f 1 1 ^ prfl ^ rj^r* a*±r f which contains accurate 
representations of the faster rhythmic activity 
uncontamina ted by the complex aperiodic processes. 

30 A polynomial is fit to the first 2m+l data points 

to determine the "trend" value at data point m+1 (i.e./ 
the middle of the range of the data points which are 
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fit with the polynomial). To determine Che* trend value 
at data point m+2 the same order polynomial is fit to 
another 2m+l data points moved one time point 
forward. This process is continued through the data 
sec until the polynomial is fit to the last 2m+l data 
points. The filter is created when the trend values 
determined by the moving polynomial are subtracted from 
the original data. Thus, the name "moving polynomial 
filter". Note that m data points at the beginning and 
end of the data set are forefeited. Since the 
procedure is equivalent to taking linear combinations 
of observations of physiological activity with 
coefficients/ the coefficients need to be calculated 
only once for a polynomial of a specified order and 
number of points and the coefficients may be used in 
all subsequent applications (see ff.G. Kendall, Time 
Series, 1973, Griffin) . 

The cubic polynomial may be expressed in the 
general form: a Q + a x t + a 2 * 2 + a 3 T 3 where t is time 
and the constants "a" are determined by the principle 
of least squares to minimize 

2 3 2 

— - -a -at-at - at) 
t=-m t 0 ~1 2~ "3 

where Y t represents the raw untransf ormed data point at 
time t, m is the number of time windows on each side of 
the midpoint of the polynomial. In our 21 point 
example, m would be ten 500 msec time windows on each 
side of the midpoint. Since we are only interested in 
generating a transformed "trend" value at time t=0 
(i.e., the midpoint of the polynomial) using the 
information of -m and +m data points, we need solve 
only for a Q which represents the weighting coefficients 
to be applied to each of the 2m+l data points to 
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produce trie best cubic polynomial fit/ i.e./ the least 
squares difference between the raw data and the 
smoothed "trend" at time t=0.. Note that a 0 is a set of 
weighting coefficients, one for each of the 2m+l data 
5 points. Through a series of mathematical calculations 
the constant a Q is calculated for a specific polynomial 
defined by a given order/ p/ and number of points/ 2m+l 
M.G. Kendall in chapter 3 of Time Series / 1973/ 
Griffin/ provides the mathematical examples for 
10 calculating the set of coefficients defining "a 0 ". The 
generalized expression to minimize the least squares 
difference for any polynomial of order p and number of 
time points 2m+l is presented below. 

l(Y - a - a t - ...a t ) 
t-m t 0 1 p 

15 The moving polynomial has two important properties. 
First/ the weights of the coefficients sum to unity. 
This is easily seen if the coefficients are applied to 
a series consisting simply of a constant which is 
repeated/ since the average must be that same 

20 constant. Second/ the weights ae symmetric about the 
middle value. 

A:s~ sta ted-above /- the— weighting. _co e.f .f ,j._ c 1 en t s used 
in the moving polynomial filter/ are those whxch 
minimize the least squares difference between the 

25 "real" observed data set/ Y^/ and the polynomial of 

order "p 1 *- and~-2nH-l— daita ipolars . In the above example / 
p=3 (i.e./ cubic) and 2ra+l=21. Once the weighting 
coefficients are calculated for a polynomial fit on any 
data set (real or hypothetical) it is not necessary to 

30 generate new coefficients. 



The following is a computer program . which 
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generates the coefficients for any polynomial of 5tft 



points. 

PROGRAM TRY ( INPUT / OUTPUT ) 
5 DIMENSION A (100) 
READ *,N1,N2 
CALL FOLYMOV (N1/N2/A) 
PRINT 50/ (A(I) /I=l/N2) 
50 FORMAT (10F8.5) 
10 STOP 
END 

SUBROUTINE POLYMOV (N1,N2,A) 

COMPUTES MOVING POLYNOMIAL FILTER COEFFICIENTS FOR C 
POLYNOMIALS OF DEGREE Nl (2 .LE. Nl .LE . 5) BASED 
15 C ON A TOTAL OF N2 POINTS. 

DIMENSION A(N2) 

IF(N1.LT.2 •OR. Nl.GT-5) RETURN 

IF(N1-KL .GT. N2) RETURN 

IF(2*tN2/2) .EQ. N2) N2=N2-1 
20 N3=N2/2 

AO=N2 

A2=0- 

A4=0 . 

A6=0- 
25 A3*0. 

DO 1 1=1, N3 

AI=I 

A2=A2+2.*AI'*2 

A4=A4+2.*AI**4 
30 A6=A6+2.*AI**6 

1 A8=A8+2.*AI**8 

IF(N1 .GT* 3) GOTO 45 

AJ=-N3-1 

DEN=AO*A4-A2*A2 
35 DEN=1./DEN 



order or lover with no limitations on the number of 
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DO 2 1=1 ,N2 
AJ=AJ+1. 
AJ2=AJ*AJ 

2 A(I)=DEN*(A4-A. J2*A2) 

5 RETURN 

45 DEN=AO*A4*A8+2.*A2*A4*A5-A6**3-AO*A6*A6-A2*A2*A8 

DEN-1 ./DEN 
C1=A4*A8-A6*A6 
C2=A4*A6-A2*A8 
10 C3=A2*A6-A4*A4 
AJ=N3-1 
DO 3 1=1 /N2 
AJ=AJ+1- 

3 A( I)=DEN*(C1+C2*AJ*AJ+C3*AJ**4) 
IS RETURN ^ 
END 

Examples are provided for It 21 and 35, point: 
polynominals of Che third order- Note chat the sum of 
Che coefficiencs approximate unity. They do not 
20 exactly sum to one since the computer approximates a 
fraction as a decimal. 

3rd order. 7 points 

-.09524 .14286 .28571 .33333 .28571 .14286 ,09524 



25 



30 



3rd order. 
-.03707 
-0.3703 
.02190 
.04717 
.06437 
.04717 
-.00442 



35 points 



-.02548 -.01460 -.00442 .00505 .01383 
.02927 .03594 .04191 

.05174 .05560 .05876 .06121 .06297 .06402 
.06402 .06297 .06121 .05876 .05560 .05174 
.04191 .03594 .02927 .02190 .01383 .00505 
-.01460 -.02548 -.03707 
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3rd order. 21 points 

-.05590 -.02484 .00294 .02746 .04871 .06669 .08140 
.09284 .10101 .10592 .10755 .10592 .10101 .09284 
.08140 .06669 .04871 .02746 .00294 -.02484 -.05590 

FIGURE 3 shows a graphical representation of the 
moving polynominal filter as it acts on the data sets 
30 output from the time window 20 of FIGURE 2a. The 
filter 90 of FIGURE 2 is constructed so as to present 
to the data 30 a set of 21 coefficients derived from a 
least squares cubic polynomial fit. The coefficients 
for each of 21 (2nH-l) points are plotted in the curve 
shown in the FIGURE 3.. Thus the actual construction of 
the curve shown in the FIGURE 3. 

Once the coefficients have been, calcuated the data 
set is multiplied by xts corresponding calculated 
coefficient to produce a template which is the smoothed 
trend as shown in FIGURE IB. The template of signals 
is then substracted from the original data set of 
FIGURE Ia which includes the data set outputs from* the 
time window in order to obtain the signal shown in 
FIGURE 1C as an output from the moving polynominal 
filter- Th_e_ value for the template associated with 
point 11 is determined by multiplying each of the 
points 1 and 21 by its corresponding normalized 
coefficient and adding the values together to obtain 
the value for the template point corresponding to the 
point 11. It is to be noted that because the 
polynominal coefficients sum to unity the normalization 
is automatically taken into account so that the value 
for the smooth template which corresponds to the raw 
data point 11 is obtained by merely adding together all 
of the multiplied values of the raw data* points 1 
through 21 with their corresponding coeffients from the 
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curve. After the summation has taken place to 
establish a corresponding template point for the raw 
data point 11 then the filter, so to speak, moves to 
data point 12 and the calculation for the data point 12 
is determined by multiplying the values of the data 
points 2 through 22 by the corresponding coefficients 
which are now centered at data point 12. Thus the 
nomenclature "moving polynomial". It can also be seen 
Chat if a 21 point polynomial filter is used then' there 
is no calculation for the points 1 through 10 nor for 
the last 10 points of the data set because the 
calculation for each particular point of raw data 
requires 10 information points from the past and 10 
information points from the future data points. The 
15 construction of a filtering device in accordance with 
the generalized polynominal equation of its 
corresponding calculated coefficients based upon the 
number of data points can be a hardware implementation 
once the values of the coefficients are- determined/ 
since they are constants for polynominals of a specific 
order and number of points. 



10 



20 



One of the key areas in thi3 polynominal filter 
which provides for a passage of useful information is 
the choosing of tFe number of- data -points which are to 

25 be fit by the polynominal. For instance in the heart 
rate embodiment if the respiration frequency which is 
the frequency of interest with regard to its effect on 
the heart rate is approximately .25 Hertz or a breath 
every 4 seconds, then in order to insure that the 

30 variance of the breathing process of interest is passed 
unattenuated the "twice duration rule" comes into 
effect. This rule basically states that in order to 
deal with both the shape of the transfer function of 
the filter which determines the shape of the frequency 
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cuc-off of the filter and the fact that the frequency 
of physiological processes are not constant but vary 
within a predictable range the use of a polynominal 
having a duration fitted locally at approximately twice 
the modal duration of the periodic process being 
studies insures that the variance of the periodic 
process of interest in is passed unattenuated or 
unembellished . Unsing this rule then the time duration 
of the polynomial filter necessary to insure the 
passage of all points of interest in a respiratory 
frequency band which is characterized by a modal 
breathing frequency of 15 breaths per minute is 8 
seconds. The breathing frequency may not be constant 
at .25 Rz but may vary from approximately .15-. 40 
Hertz. The modal frequency (i.e./ dominant peridocity) 
may vaxy among subjects/ but most adults breathe with a 
modal frequency between .15 and .40. Therefore/ it can 
be seen that a 21 point polynomial filter generally 
satisfies this requirement as each point taken from the 
time window is separated by 500 msec and therefore the 
21 points encompass a time frame of 10.5 seconds which 
is sufficient to pass those events of interest 
occurring because of the respiration. It should also 
be noted that in some instances depending upon the 
application- of the -filtering device it is not 
advantageous to increase the number of points too 
broadly because this would effectively pass the 
variance associated with slower periodic and aperiodic 
processes. Since physiuological processes are not 
perfect sine waves, passing lower frequency 
physiological activity will always result in the 
passing of higher frequency harmonics which may be 
superimposed on the variance of the frequency band of 
interest. 
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She selection of Che cubic order for the 
polynomial considerably enhances Che ability of the 
filter to respond to and "fit" cubic trends. Higher 
order polynomials may, on the local level/ start to map 
5 into the faster periodic activity/ while lower order 

polynomials may not *bend n with the baseline trend. In 
summation then, the cubic order polynomial is chosen 
because it maps the dips in the response pattern 
without mapping into tne faster periodic component. 

10 a bandpass filter 94 is added to the output of the 

MPP in order to reduce the outpuc variance to only 
periodicities of interest. This necessitates prior 
knowledge of the physiological rhythms which are bing 
studied. For instance, with regard to respiratory 

15 sinus arrhythmia (RSA), it is determined by the natural 
range of breathing for the age for the human being 
tested. The output of the bandpass filter is 
equivalent to the variance of the frequency band of 
interest and is output to a variance calculator and 

20 Display 76. Because all physiological oscillations are 
not perfect sine waves/ the periodicities are 
manifested across a band of frequencies. Thus, the sum 
of the variances associated with this band of 
.frequencies is desired. This sum may be calculated per 

25 time period with spectral analysis a^id by summing- t he- 
spectral densities associated with the frequencies of 
interest. Likewise/ in situations in which rapid and 
continuous estimates are required, the. variance may be 
calculated witn. traditional descriptive statistics from 

30 the output of the bandpass filter. When the rhythmic 
oscillation in heart period is clearly distinguished 
from the background activity through the above 
processes/ the variance which is calculated by either 
of the above two methods provides an estimate of the 
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vagal cone* 

In the above analysis the bandpass filter 94 which 
takes the output of the MPF passes only the frequencies 
of interest which in the case of RSA with adults 
requires a bandpass of 0.15-0,40 cycles per second in 
order to allow only the periodicities associated with 
breathing to pass. 

When the coupling of two or more physiological 
processes needs to be evaluated, the simultaneous 
output of two series fed through the filtering 
procedure may be analyzed with cross-spectral analysis 
and a weighted coherence may be calculated ♦ This is 
shown in FIGURE 2 as block 98. Since the weighted 
coherence (see Porges et al., "New time-series 
statistic for detecting rhythraxc co-occurrence in the 
frequency domain: The weighted coherence and its 
application to psychophysiological research, " 
Psychological Bulletin , 1980/ 88/ 580-587) is derived 
by weighting the coherence function across a band of 
frequencies by the spectral densities/ accurate 
estimates of the spectral densities provided by the 
above procedures are necessary to generate an 
appropriate -estimate of the weighted coherence. 

The FIGURES 4 a,b shows a clinical application of . 
the results of the analysis detailed above by the 
structure of FIGURE 2. The detection and evaluation of 
the component of the heart rate pattern which is 
directly influenced by the central nervous system/ 
called the vagal tone, is detected by the output of the 
FIGURE 2. Aji evaluation of the vagal tone in a group 
of normal infants and a group containing a variety of 
clinical pathologies including severe brain damage is 
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detailed in the FIGURE 4b. The total heart period 
variability was also collected. There is an apparent 
monotonia relationship between the vagal tone detected 
by the embodiment employing the structure of FIGURE 3 
and the severity of clinical dysfunction- Then the 
same infants were ranked in terms of their heart period 
variability as in FIGURE 4a/ there was a clear 
distinction only between those who died and all other 
infants. While heart period variability clearly 
distinguished between brain death infants/ with their 
characteristic absence of neural influence on the 
heart/ and all other infants/ it did not distinguish 
among the various infants having neural tube defects, 
respiratory distress syndrome/ and normal infants. The 
letters on the Scale of FIGURES 4a /b represent the 
diagnosis or insult associated with individual infants 
(A » asphyxia; BPD = bronchio-pulmonary dysplasia; CA » 
cardiac arrest; H = hydrocephalic r M - microcephalic; 
N = normal; RDS = respiratory-distress syndrome; 
* = denotes infants who subsequently expired). Heat 
period did not reliably discriminate among the various 
pathologies although there was a tendency among the 
severly brain damaged to have ajiort heart peiods. 
Categorization with the heart period varia&jrlj^ ty o gX 
FIGURE 4a partitioned the infants - into two global 
categories while classification by vagal tone in FIGURE 
4b as measured with the FIGURE 2 embodiment allowed for 
a continuum, of severity of neuopathology . Thus/ 
although the heart period variability of FIGURE 4a, 
which was used previously as an indicator in this type 
of diagnosis/ is sensitive to gross dysfunction/ it is 
quite clear that the vagal tone used in FIGURE 4b is a 
more sensitive index to individual differences in 
central dysfunction . 
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Obviously/ numerous modif ication. and variations of 
the present invention are possible in light of the 
above teachings. It is therefore to be understood that 
within the scope of the appended claims/ the invention 
may be practiced otherwise than as specifically 
described herein. 
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Claims 

1. A system for detecting amplitude variations in 
the rhythmic oscillations of heart period in a 
predetermined frequency band of interest, comprising: 

sensor means for detecting an occurrence of an 
event in the cycle of a heartbeat; 

timing means for determining intervals between 
said event in successive heartbeats; 

buffer means for storing said intervals between 
successive beats; 

means- for r ea ding out values stored in said buffer 
means at predetermined intervals; 

first filter means responsive to the output of 
said means for reading out values including a means for 
determining an aperiodic portion of said read out 
signal from said buffer means and subtraction means for 
subtracting said aperiodic portion of said read out 
signal from said read out signal of said buffer means; 

second" filter- means -for- receiving, the piJtp_ut of 
said first filter means and outputting a signal in a 
band-pass region which band— pass region is determined 
by said predetermined frequency band of interest; and 

output means for calculating the variance of said 
rhythmic oscillation. 

2- A system for detecting amplitude variations in 
the rhythmic oscillations of a physiological response 
pattern occurring in a predetermined frequency band of 
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interest/ comprising: 

sensor means for detecting said physiological 
response pattern: : 

analog-digital converting means for converting 
said response pattern into digital form whereby the 
sampling frequency of said analog to digital converter 
is determined by said frequency band of interest; 

first filter means responsive to the output of 
said analog to digital converter including an aperiodic 
determining means for determining the aperiodic portion 
of said physiological response pattern from said analog 
to digital converter and subtraction means for 

subtracting said aperiodic portion of said output of \ 
said analog-digital converter from said output of said 
analog— digital converter ; 

second filter means for receiving the output of 
said fxrst filter means for outputting a signal in a 
band-pass region which band-pass region is determined 
by said predetermined frequency band of interest; and 

output means for calculating the variance of said 
rhythmic oscillation. 

3. A method- for. detecting amplitude variations in 
the rhythmic oscillations of heart period in a 
predetermined frequency band of interest comprising the 
steps of: 

detecting the occurrence of an event in a cycle of 
a heartbeat; 
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deterraining Che incervals between said event in 
successive heartbeats; 

storing said determined intervals; 

reading out said stored values at predetermined 
5 intervals; 

performing a first filtering function on said 
readout values including removing the aperiodic portion 
of said readout signal; 

receiving the physiological response pattern 
10 signal with its aperiodic portion removed and 

outputting a signal in a band pass region which band 
pass region is determined by said predetermined 
frequency band of interest r and 

calculating the variance of said rhythmic 
15 oscillation . 

4. A method for detecting amplitude variations in 
the rhythmic oscillations of a physiological response 
pattern occurring in a predetermined frequency band of 
interest/ compj:isi n 9 the steps of: 

2o detecting said physiological response pattern; 

converting said response pattern into digital form 
by sampling said response pattern at a frequency 
determined by said predetermined frequency band; 



25 



performing a first filtering function on said 
digital form of said response pattern including 
removing the aperiodic portion of said digital form of 
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said response pattern; 

receiving said digital form of said response 
pattern with its aperiodic portion removed and 
outputting a signal in a band— pass region which band-^ 
pass region is determined by said predetermined 
frequency band; and 

calculating the variance of said rhythmic 
oscillation. 

5. The system according to Claim 1 wherein said 
predetermined frequency band of interest is associated 
with respiratory activity. 

6. The system according to Claim 1 further 
comprising a compensator means for receiving the output 
of said buffer means and adjusting said intervals 
between said event in successive heart beats to within 
a predetermined range. 

7. A system according to Claim 1 wherein said 
first filter means comprises a moving polynomial, filter 
constructed on the basis of a 2ra+l point polynomial of 
order n whereby said aperiodic portion of said readout 
signal is determined and subtracted from said readout 
signal to produce a residual signal to be outputted to 
said second filter means. 

8. The system according to Claim 7 wherein m-10 
and n=3. 

9. The system according to Claim 7 wherein m=25 
and n=3. 
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10. The system according no Claim 2 wherein a 
plurality of physiological response patterns are 
detected and wherein said output means further includes 
a coherence measurement means for providing a weighted 

5 coherence between the calculated variance of each of 
said response pattern occurring in said predetermined 
frequency range of interest. 

11. The system according to Claim 2 wherein said 
predetermined frequency band of interest is associated 

10 with respiratory activity. 

12. The system according to Claim 2 further 
comprising a compensator means for receiving the output 
of said buffer means and adjusting said intervals 
between said event in successive heart beats to within 

15 a predetermined range. 

13. A system according to Claim 2 wherein said 
first filter means comprises a moving polynomial filter 
constructed on the basis of a 2m+l point polynomial of 
order n whereby said aperiodic portion of said readout 

20 signal is determined and subtracted from said readout 
signal to produce a residual signal to be outputted to 
said second filter means. 

14. The system according to Claim 13 wherein m=10 
and n=3 - 

15. The system according to Claim' 13 wherein m=25 
and n=3 . 

16. The method according to Claim 3 wherein the 
step of storing said determined intervals includes the 
seep of compensating for said determined intervals 
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which are outside of a predetermined range so that all 
intervals are within said predetermined range. 

17. The method according to Claim 3 wherein the 
step of performing a first filtering function includes 

5 the step of determining the periodic portion of said 
readout signal through the use of a moving polynomial 
filter based upon a polynomial equation having 2ra+l 
point and being of the nth order. 

18. The method according to Claim 17 where m=10 
10 and n=3. 

19. The method according to Claim 17 where m=25 
and n-3 . 

20. The method according to Claim 4 wherein the 
step of storing said determined intervals includes the 

IS step of compensating for said determined intervals 

which are outside of a predetermined range so tha.t all 
intervals are within said predetermined range. 



20 



21. The method according to Claim 4 wherein the 
step of performing a first filtering function includes 
the step of determining the periodic portion of said 
readout signal through the use of a movirrtg- polynomial 
filter based upon a polynomial equation having 2m+i 
points and being of the nth order. 

22. The method according to Claim 21 where m=10 
25 and n=3. 

23. The method according to Claim 21 where m=25 
and n=3. 
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24, A system for detecting amplitude variations 
in rhythmic oscillations of a point process in a 
predetermined frequency band of interest/ comprising: 

sensor means for detecting an occurrence of a 
5 plurality of events in the cycle of each of said, 
rhythmic oscillations ; 

timing means for determining intervals between the 
occurrence of successive of said events; 

buffer means for storing said intervals; 

10 means for reading out values stored in said buffer 

means at predetermined intervals; 

first filter means responsive to the output of 
said means for reading out values including a means for 
determining an aperiodic portion of said read out 
15 signal from said buffer means and subtraction means for 
subtracting said aperiodic portion of said read out 
signal from said read out signal of said buffer means; 

second filter means for receiving the output of 
said first filter means and outputting a signal in a 
20 band-pass region which band-pass- region as determined 
by said predetermined frequency band; and 

output- means f or~ calculating the variance of said 
rhythmic oscillation. 

25. A system for detecting amplitude variations 
25 in the rhythmic oscillations of a continuous process in 

a predetermined frequency band of interest/ comprising: 
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sensor means for detecting and evaluating said 
continuous process and oucputting a signal; 

analog to digital converting means for converting 
said output signal of said sensor means into digital 
form whereby the sampling frequency of said analog to 
digztal converter is determined by said frequency band 
of interest; 

first filter means responsive to the output of 
said analog to digital converter including an aperiodic 
portion determining means for determining the aperiodic 
portion of said output of said analog to digital 
converter and a subtraction means for subtracting said 
aperiodic portion of said output of said analog to 
digital converter from said output of said analog to 
digital converter; 

second filter means for receiving the output of 
said first filter means for outputting a signal in a 
band-pass region which band-pass region is determined 
by said predetermined frequency band of interest; 

output means for calculating the variance of said 
rhythmic oscillation * 

26* A method for detecting amplitude variations 
in the rhythmic oscillations of a continuous process in 
a predetermined frequency band of inceceat^compcising 
the steps of: 

detecting said continuous process ; 

converting said detected process into digital form 
by sampling said process at a frequency determined by 
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said predetermined frequency band; 

performing a first filtering function on said 
digital form of said process including removing M 
aperiodic portion of said digital form of said process; 

5 receiving said digital form of said process with 

its aperiodic portion removed and outputting a signal 
in a band-pass region which band-pass region is - 
determined by said predetermined frequency band; and 

calculating the variance of said rhythmic 
10 oscillation. 

27. A method- for detecting amplitude variations 
in the rhythmi.c oscillations of a point process in a 
predetermined frequency band of interest comprising the 
steps of: 

ls detecting the occurrence of an event in said point 

process ; 

determining the intervals between successive 
events in said point process; 

storing said determined intervals: 

2Q reading out said stored values at predetermined 

i nutfi. va-l sr 

performing a first filtering function on said read 
out values including removing the aperiodic portion of 
said read out signal; 



25 receiving the point process signal with its 
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aperiodic portion removed and outputting a signal in a 
band-pass region which band— pass region is determined 
by said predetermined frequency band of interest; 

calculating the variance of said rhythmic 
oscillation. 

28. A system for detecting amplitude variations 
in rhythmic oscillations of a plurality. of signals in a 
predetermined frequency band wherein each of said 
rhythmic oscillations occurs in one of a point process 
and a continuous process T comprising: 

first sensor means for detecting an occurrence of 
a plurality of events in the cycle of each of said 
rhythmic oscillations occurring in a point process; 

second sensor means for detecting and evaluating 
said continuous process of each of said rhythmic 
oscillations occurring in a continuous process; 

timing means responsive to the output of said 
first sensor means for determining intervals between 
the occurrence of successive of said events; 

buffer means for storing said intervals; 

means for reading, out- said ^f^rv ala stored in 
said buffer means at predetermined intervals; 

analog to digital converting means for converting 
the output of said second sensor means into digital 
form whereby the sampling frequency of said analog to 
digital converter is determined by said frequency band 
of interest; 
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first: filter means responsive co the output of 
said analog to digital converter means and said means 
for reading out values including a means for 
determining an aperiodic portion of at least one of 
said read out intervals and the output of said analog 
to digital converter means and subtraction means for 
subtracting said periodic portion of at least one of 
said read out intervals and said output of said analog 
to digital converter from at least one of said read out 
intervals of said buffer means and said output of said 
analog to digital converter; 

second filter means for receiving the output of 
said first filter means for outputting a signal in a 
band-pass region which band-pass region is determined 
by said predetermined frequency of interest; 

output means for calculating the variance of each 
of said rhythmic oscillations; and 

a coherence measurement means for providing a 
weighted coherence between the calculated variance of 
any two of said plurality of rhythmic oscillations- 
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